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Andreev levels deplete energy states above the superconductive gap, which leads to the peculiar non- 
monotonous crossover in the local density of states of mesoscopic superconductor/normal-metal/superconductor 
junctions. This effect is especially pronounced in the case when the normal metal bridge length L is small com- 
pared to the superconductive coherence length £. Remarkable property of the crossover function is that it 
vanishes not only at the proximity induced gap 6^ but also at the superconductive gap A. Analytical expressions 
for the density of states at the both gap edges, as well as general structure of the crossover are discussed. 

PACS numbers: 74.45. +c 
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Experimental advances in probing systems at the meso- 
scopic scale revived interest to the proximity related 
problems in superconductor - normal metal (SN) heterostruc- 
tures.* The most simple physical quantity reflecting proxim- 
ity effect is the local density of states (LDOS) p{e, r), which 
can be measured in any spatial point r at given energy e us- 
ing scanning tunneling microscopy. The effects of supercon- 
ductive correlations on the spectrum of a normal metal are 
especially dramatic in restricted geometries. For example, 
in the case of superconductor-normal metal-superconductor 
(SNS) junction, proximity effect induces an energy gap eg 
in excitation spectrum of a normal metal with the square 

-ije/eg - 1 in the density of states 
<K €„ Ref. 



root singularity p(e, r) 
just above the threshold e - £g ^ eg 
(here and in what follows, p will be measured in units of the 
bare normal metal density of states v at Fermi energy). The 
most recent theoretical interest was devoted either to meso- 



scopio 



13.14.16.17 



or quantum 



18.19.20.21 



fluctuation effects on top 



of mean-field result a^i^i^''°i^'''^ that smear hard gap below £g 
and lead to the so called subgap tail states with nonvanishing 
p oc exp [- g(l - e/et,)*^"'''^"*] at - e < eg, where g is the 
dimensionless normal wire conductance and d is the effective 
system dimensionality. The latter is essentially a nonpertur- 
bative result that requires instantonlike approach within cr- 
mode l'^'^" or relies on methods of random matrix theory ji^ii^ 
Surprisingly, after all of these advances, there is something 
interesting to discuss about proximity induced properties of 
the SNS junctions even at the level of quasiclassical approxi- 
mation by employing Usadel equations.-^ The purpose of this 
work is to point out a subtle feature of the crossover in the lo- 
cal density of states of mesoscopic SNS junctions. The latter 
was seen in some early and recent studieS j^'"''^''" ^''^ however, 
neither emphasized nor theoretically addressed. 

To this end, consider normal wire (N) of length L and 
width W located between two superconductive electrodes 
(S). In what follows, we concentrate on diffusive quasi-one- 
dimensional geometry and the short wire limit L <k ^, where 
^ - y/Ds /A is superconductive coherence length, with Ds as 
the diffusion coefficient in the superconductor and A as the en- 
ergy gap (hereafter, ^ = 1). The center of the wire is chosen to 
be at jc = and boundaries with superconductors at x = +L/2 
correspondingly, where x is the coordinate along the wire. 

Under the condition L «; ^, the proximity effect is es- 




FIG. 1 : (Color online) Schematic plot for the local density of states 
crossover in the short L <c ^ diffusive SNS junction. 



pecially strong — superconductive correlations penetrate the 
entire volume of the normal region. As the result, the in- 
duced energy gap in normal wire eg is large and turns out to 
be of the same order as the gap in the superconductor itself 
eg - A - 6 A, with the finite size correction 6 A ~ A^ /e^j^ <K A, 
where ejh = D^IL^ is the Thouless energy and is the dif- 
fusion coefficient in the normal bridge. The latter should be 
contrasted to the long wire limit L » ^, where the proxim- 
ity effect is weak and the induced gap is eg ~ ejh <s: A. Just 



above the gap e - « eg, density of states has square root 
singularity, which is similar to the L » ^ case 



7.8.9.10.11.12 



see 



Fig. [T] which is a robust property of quasiclassical approx- 
imation. However, the prefactor for L <k ^ is significantly 
enhanced p(e, x) oc (ej-/,/A)^ -y/e/e^ - 1 (note here that propor- 
tionality sign implies characteristic energy dependence, while 
the exact numerical coefficient is different for a given coordi- 
nate X along the wire). The value of eg is a property of the 
spectrum, thus it is x independent. 

One would naturally expect that above the proximity in- 
duced gap eg, local density of states goes through the max- 
imum pmax - p(e - A, x) and then crosses over to the BCS 



2 



like DOS pbcs - e/ Ve^ - at e > A, which finally satu- 
rates to unity p — > 1 at e — > oo. Surprisingly, however, the 
crossover scenario is different. The density of states indeed 
reaches the maximum, which occurs at e ~ A - (5A/2, but then 
decreases and vanishes to zero at superconductive gap A with 
quarctic power law behavior p(e, x) oc {eTi,/A)^^^-\/\e/A - 1| 
for |e - A| < 6A. Finally, p(e, x) grows back for e > A and at 
the energy scale e* ~ A + 6A crosses over to the pscs ■ Indeed, 
observe that pscsie*) ~ er/i/A ~ p{e*,x), see Fig. [T] It is 
important to mention here that the discussed feature appears 
only at the level of local density of states. Energy dependence 
of the global density of states, which is integrated over the 
entire volume, does not show dip at e = A (this point will be 
discussed later in the text). 

It seems that crossover picture presented in the Fig. [1] was 
already numerically seen in previous studies, "''^''"^ '^ how- 
ever the origin of the soft gap at e = A was never addressed. 
One should also note that numerics was performed always for 
the not too short junctions L > ^. It is certainly unusual to 
find zero in the density of states at the gap edge A. In order to 
get some insight into this bizarre feature, let us consider a toy 
model, which crudely mimics the system under consideration 
and gives some hints for qualitative understanding. Quantita- 
tive theory of the outlined crossover will be developed follow- 
ing this discussion. 

Imagine chaotic quantum dot (QD) sandwiched between 
two superconductors. A quantum dot is really a zero dimen- 
sional system, in the sense that L/^ — > (or equivalently, 
Alejh — > 0). However, one will momentarily see that finite 
tunneling rate into the dot F plays an effective role of Thouless 
energy ejh- By following Ref. 23 one may construct scatter- 
ing states and determine quasiparticles excitation spectrum. 
If in addition we assume that QD supports only one trans- 
verse propagating mode, then the discrete spectrum obtained 
from the poles of the scattering matrix Sie) consists of a single 
nondegenerate Andreev state at energy e [0, A], satisfying 
Q.{e) + Fe^ VA^ - - , where the function Q(e) is defined 
by Q(e) = (A^ - e^)(e^ - T^/A-). The density of states in the 
superconductor-quantum dot-superconductor system is given 
by p(e) = pBcs(e) + 5p{e), where 



6p(e) 



1 d 
2ni de 



In 



'Q(e) + iTe^ ^Je^ - A^] 



,Q(e) - /re2 Ve2 - A^ ) 



(1) 



which follows from the relation p(e) = Pbcs + 
(1 /27r0(d/de) lnDet.S(e) between DOS and the scattering ma- 
trix.— It is easy to check that in the limit F » A Andreev 
level is positioned at eo A - SA-'/F^, which resembles the 
expression for eg in the case of the diffusive wire discussed 
above, where F indeed plays the role of Thouless energy. The 
presence of an Andreev level changes density of states p(e) 
in two ways. The first contribution 5p\{e) originating from 
d In Det>S/de term of Eq. ([1]) belongs to the subgap part of the 
spectrum e € [0, A] and has structure of the form 



^ ^ ^ Id e-eo-iO 
Spi(e) = z— ;5-ln — 



^S(e-e„), 



(2) 



which is nothing else but DOS associated with the single An- 
dreev state. Most interestingly, Andreev level changes p(e) 



above the superconducting gap as well. Indeed, it follows 
from the Eq. ([T]i that at energies e > A the density of states 
p(e) gets the correction 



Sp2ie) = -PBCsie)- 



F (2A2 - e2)r2/4 - e^{e^ + A^) 



^ (e2 - A2)(e2 - ry4y + FV 



(3) 



Observe that in the immediate vicinity of the superconductive 
gap e - A< A^/T^ the coiTection 6p2{e) ^ -(r/4n:A^)pBcs(e) 
is negative, implying that Andreev level suppresses bulk su- 
perconductive density of states pscs (f) above the gap. At en- 
ergy e* ~ A + A^/F^ the correction given by Eq. (O reaches 
its maximum, while pscs is recovered at large energies e > F, 
where 6p2 decays as 6p2{e) ~ F/ne^. Based on this example, 
it appears that Andreev level tends to deplete bulk BCS den- 
sity of states at energies e - A < A^/F^. One should note that 
the Andreev level leads to pure redistribution of the energy 
states — there are no additional states above the gap; indeed, 
^P2(f)de = 0. If QD supports not only one but a large 
number of the transverse propagating channels, then cummu- 
lative negative dpiie) may compensate pscsi^), which leads 
to the vanishing density of states pie) at the gap edge e - A, 
see Fig.[T] 

Having discussed the qualitative picture of the crossover, 
let us now turn to the quantitative description. The quasi- 
classical approach to diffusive SNS structures is based on the 
Usadel equation^^ for the retarded Green's function S^(r, e). 
For the latter, we will employ angular parametrization^^ = 
T, COS 6 + T-c sin cos </> + t,j sin 9 sin (f>, where r,- is the set of 
Pauli matrices. In the absence of the phase difference between 
the S terminals, one can set = and Usadel equations for 
quasi-one-dimension SNS geometry acquires the form 

'DN{(f6N/dx^) + liesinON = |x| < L/2 

Ds {d^6s /dx^) + lie sin Os + 2 A cos 6s ^ \x\ > L/2 ^ ' 

where 0n{S){£, x) are the Green's function angles in N(S) parts 
of the junction, correspondingly. In writing Eq. (|4]i, we as- 
sumed step function pair-potential A(x) = Ari{\x\ - L/2), with 
ri(x) = 1 if jc > and ri(x) - otherwise. The applicability of 
this approximation relies on the condition that the width W of 
the junction is small compared to the coherence length. In this 
case, nonuniformities in A(r) extend only over the distance of 
order W from the junction, which is due to the geometrically 
constrained influence of the narrow junction on the bulk su- 
perconductor 

In the absence of additional tunnel barriers at SN- 
interfaces, Eqs. (|4^) and ^p) are supplemented by the fol- 
lowing boundary conditionsi^ 



0iv(e, x)L-=±z,/2 = 0s(e,x)l^ 



=±£/2 : 



ctn- 



dOi^ie, x) 



dx 



x=±L/2 



as- 



d0s(e,x) 



dx 



x=±L/2 



(5a) 



(5b) 



where ctn^s) are the conductivities of normal metal and super- 
conductor Knowing solutions of Usadel equations, one finds 
local density of states from the general expression 



p(e, x) - Re[cos 6{e, x)]. 



(6) 
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For future convenience, we rotate the Green's function an- 
gles as Of^ - n/2 + /j?^ and 0s = Obcs + i^s, where Obcs - 
7t/2 + /arsinhy with 7 = e/ Vl - e^, and introduce dimension- 
less variables e = e/A, A - x/L. After the rotation, Usadel 
Eq. (|4j)) for superconducting sides of the junction becomes 
real and can be easily integrated, providing 



§s(£,'i) = 4artanh{exp[-(|i| - l/2)/^s]] 



(7) 



where we have introduced the energy depending supercon- 
ductive coherence length = "sI ItiJ)" ~ Equation 
(|4^) may also be exactly solved in terms of elliptic functions; 
however, this exact solution is not needed for our purposes. 
Indeed, observe that it follows from the boundary condition 
[Eq. (|5^)] that in the energy range 1 - s 1, the normal 

metal phase Re[i?iv(e, /I)] oc arcsinhy x In -yj-^ » 1 is large 
everywhere in the N part of the junction. Thus, one may ap- 
proximate sin^Af =s exp(??Af)/2 and solve Eq. ^) in closed 
form 



&N(e, /I) = &o(b) - In [ cosh2(^/^^)] . 



with normal side coherence length being defined as - 
yj^^ and )?o = '?A'(e, 0). We now introduce uq - expi&o and 

Us = exp(arcsinh7 + §^), where i?| = §s (s, ± 1 /2), to rewrite 
the boundary condition [Eq. (jSJ))] as 



usly-2 = xy(uQ - Us) , 



(9) 



where the interface parameter x - cr^Ds /ct^Dn measures 
the mismatch of conductivities and diffusion coefficients at 
the SN boundaries. By using solutions ^ and (O, together 
with Eq. (|9]l, one eliminates the unknown uq and arrives to the 
algebraic equation for z - us ly — 2, 



yeA 



(10) 



where the single parameter k - xy^ scaling function 'Jiz, k) is 
given by 



z + (z + 2)k 



arctanh 



z + {z + 2)k 



(11) 



Knowing the solution of Eq. fTOl l, one finds the density of 
states at the SN interfaces as 



p(6, ±1/2)= I Im[z(e)]. 



(12) 



At the same time, us{e) together with Eqs. ©-([H]) provide an 
explicit information about the local density of states p(e, x) at 
any position x along the wire. 

By looking at Eq. (ITOb . one sees that its right hand side 
grows to infinity when e — » 1, while its left hand side has 
an absolute maximum for certain value of z- This implies 
that for all energies below some threshold e^, Eq. ([TOl i has the 
only real solution for z providing pie) = as it follows from 



Eq. (ITST i. The condition e - Eg when Eq. (ITOb has a complex 
solution for z for the first time defines the proximity induced 
energy gap Eg. For energies above the gap e > Eg, the density 
of states is nonzero since Im[z] 0. It turns out that Eq. (ITOl l 
possesses two qualitatively different solutions depending on 
the value of the interface parameter x. 

The limit of strong superconductor x <K A^/e^^^. In this 
case, /f <K 1 and 7 function determined by Eq. ( ITTT i has 
the following asymptotes: 3'(z, k) » y[zJ\K for z <K 1 and 
'S{z,K) ~ V^74iln(2/A-) for z » 1. It means that the only 
relevant z, which determine the maximum of 3^ are those 
z ~ K. In this region, 'S(z,k) may be approximated as 



nZ,K) 



arctanh 



As a result, the absolute max- 



imum 'Sm - 9^(z = Zm) occurs at point z,„ ■x, 4.5*', which cor- 
responds to 0.5. Then, the gap determining condition 
gives 



T - 



ygSgA 



E„ = I - 
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(13) 



(8) where we have used the notation yg - yisg). Just above the 



gap E-Eg «i Eg, one can expand 3^ in Taylor series around the 
maximum J k J,„ + b(z - z,„)^ with b = (l/2)idP-:f/dz\=,,„ 

to find z ~ z,„ + i ^ using now def- 

inition ([T2I 1. one finds for the density of states just above the 
proximity induced gap at SN-interface, 



p(,,,l/2)oc(^)^^ 



e - e„ <s fa , 



(14) 



where the numerical coefficient of the order of unity was omit- 
ted. For the other limiting case, in the vicinity of the super- 
conductive gap e ~ 1, Eq. ( ITOl i is solved by z ^ -n^p^{\-Aip) 

with P - y/l - E, which gives for the density of states. 



p(.,±l/2)oc(_J ^_ 



|e - A| < dA. 



(15) 



This asymptotic result holds above A as well. Observe that at 
e ~ €g+6A/2, Eqs. ST4\ and ( fTSI ) crossover to each other, while 
at e ~ A 4- 6A, Eq. dTSl l crossovers to the BCS like density of 
states. 

The limit of weak superconductor x » A^/e^^. This 
limiting case corresponds to the situation when a- » 1 and 

the expression for 5" greatly simplifies S'iz, k) ~ -^j^T^ 



At 



2, the function 3^ has the maximum 3^,„ 



VT78^ 



were 



Kg = xyg, so that the gap determining condition is different 
from Eq. ( [T3l ) and reads 




2 \eTh 



2/3 



(16) 



To calculate the asymptotes for density of states at both gap 
edges, one follows the same steps as in the previous case and 



4 



finds 



^^^'"^/2)oc(^) It^, e-e,«e„ (17a) 



p(e,±l/2)(x y^^ll^, |e-A|<M*, (17b) 

where 6A* ~ A{xA/£Thf'^- Equations ^ and ^ 
complement our qualitative considerations presented at the 
beginning of this paper 

At this point, let us discuss the obtained results and limits 
of their applicability, (i) We have studied the energy depen- 
dence of the local density of states for short (L <K ^) dif- 
fusive SNS junctions. Although p(e, x) was analytically cal- 
culated at SN interfaces only it turns out that its energy de- 
pendence is generic for any x e [-L/2,L/2] and given by 
Eqs. (O, (fTSl l and ( fTTb . The exact numerical prefactor, how- 
ever, is X dependent and should be determined numerically, 
(ii) Let us stress that the discussed feature in the p(e,x) at 
e ~ A disappears at the level of global (p(e)), which is inte- 
grated over the volume, density of states. Indeed, the spatial 



integration brings an additional factor of f^, where d is an 
eff'ective dimensionality of superconductor, which is due to 
the long spatial tails of Andreev states penetrating deep inside 
the superconductor [note that because of these tails, it is not 
correct to integrate p{e,x) over x e [-L/2, L/2] only]. For 
quasi-one-dimensional geometry discussed here, a factor of 
ij^y^'^ gained after x integration exactly compensates 
the dip in p{e,x) at e ~ A [see Eqs. ( fTSl ) and (fT7b)1. leading 
to the finite value {p(e-A)} ~ exh/A. For d > 1, the density 
of states at superconductive gap edge e ~ A should diverge as 
a certain power-law (p(e)) ~ (e - A)"P for p > 0. (iii) It fol- 
lows from the numerical analysis'' that presence of additional 
tunnel barriers at SN interfaces alters soft gap and leads to the 
nonzero density of states at the gap edge A. (iv) Finite super- 
conductive phase <p imposed across the junction shifts position 
of the proximity induced gap'' eg{<p), such that eg{<p - n) = 0, 
and also changes shape of the crossover function. However, 
zero in the LDOS at e = A and asymptotes at both gap edges 
persist. 
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